Optical Interferometers with Reduced Sensitivity to Thermal Noise 
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A fundamental limit to the sensitivity of optical interferometry is thermal noise that drives fluc- 
tuations in the positions of the surfaces of the interferometer's mirrors, and thereby in the phase of 
the intracavity field. Schemes for reducing this thermally driven phase noise are presented in which 
phase shifts from concomitant strains at the surface and in the bulk of the substrate compensate 
the phase shift due to the displacement of the surface. Although the position of the physical surface 
fluctuates, the optical phase upon reflection can have reduced sensitivity to this motion. 



Thermal noise presents a fundamental limit to mea- 
surement sensitivity in diverse areas of science and tech- 
nology pp. One important setting is that of optical in- 
terferometry in which otherwise stable structures expe- 
rience small, thermally driven fluctuations in their di- 
mensions in applications ranging from frequency metrol- 
ogy 02 S], to gravitational wave detection [5j [6], to 
the realization of quantum behavior for macroscopic ob- 
jects pp. 

The dominant limitation to length stability for some 
interferometers originates from thermally driven dis- 
placement noise for the reflective surfaces of the mirror 
substrates, and not from the external, support struc- 
ture [3j [U EJ. The fluctuations of the mirror surfaces 
are of fundamental origin and arise from dissipation 
in the elasticity of the substrate as demanded by the 
Fluctuation-Dissipation Theorem (FDT) [U [10] . 

Beyond displacements driven by thermal noise in the 
substrate itself [HJ E2 US US US US], diverse other 
sources of mechanical noise have been identified in recent 
years, including noise from frictional losses in the mate- 
rials that form the mirror coating p~4j [17] , thermoelastic- 
damping in the substrate and coating p~6j [18] [19], and 
thermo-refractive noise [20]. Measurements of phase 
noise in rigid [2] [3] [4] and suspended interferome- 
ters [HI [22] [23] have confirmed many characteristics of 
individual noise sources. 

Various avenues have been followed for reducing ther- 
mal noise in optical interferometers, the most significant 
being lowering mechanical losses for the substrate [TJ [24] 
and, more recently, the coating [23j [25]. New designs 
for advanced interferometers include the use of a corner 
reflector [26 or a short Fabry- Perot (FP) cavity [27] to 
replace the usual single surface of a mirror. 

In this Letter, we propose a new strategy for reduc- 
ing thermally driven phase noise in optical interferome- 
ters. Fundamental to our proposal is the observation that 
stochastic displacements 5u z perpendicular to the mir- 
ror's surface are necessarily accompanied by correlated 
strains in the underlying materials of the mirror coating 
and substrate. In a conventional setting, the phase shift 
5/3 due to these strains in the coating and substrate are 
small compared to the phase shift 50 from the surface 
motion 5u z . However, by suitable design, it should be 
possible to achieve a total phase shift 5$ = 5/3 + 50 ~ 




FIG. 1: (color online) (a) The undistorted shape of a sap- 
phire substrate of radius a = 1.5 mm and thickness / = 1 
mm is depicted by the shaded region and wire frame. Ax- 
ial displacements u z (x, y, z) and strains e zz (x, y, z) are shown 
across planes at z — 0,//2,/ (top to bottom) for the eigen- 
mode £ ( r ) with eigenfrequency ujq/2ti = 2.22 MHz. Plotted 
are contours for u z (x, y, z) on which e zz (x, y, z) is color coded 
at a time of maximum axial extension, where u z (x,y, z) = 
absent excitation. The phase for propagation from z = zr to 
the substrate and back is modified by surface motion u z at 
z = 0. (b) Coating stack for a high reflectivity mirror with 
embedded FP cavity for high strain sensitivity. 



for the reflected field. That is, although the physical 
surface of the mirror is subject to random displacements 
5u zi 5$ can be insensitive to these displacements with 
50 dynamically compensated by 5(3 from the coating and 
substrate. Although our current analysis is for "Brown- 
ian noise" of the substrate, our methodology should be 
applicable to thermoelastic and coating noise as well. 

As a first example, consider thermal noise for an eigen- 
mode of a cylindrical mirror with mechanical resonance 
frequency c^o- For frequencies cu ~ cuq and quality fac- 
tor Qo ^> 1, microscopic thermal noise excites the entire 
"shape" of the relevant eigenmode ^ (r), with a single, 
overall amplitude set by equipartition of energy [HI US • 
The amplitude 5u z and accompanying strain e zz for ther- 
mally driven motion perpendicular to the mirror surface 
at z = are then determined directly from £o( r )- 

Figure w[ a) illustrates a particular axisymmetric eigen- 
mode ^ (r ) for a sapphire substrate of mass M as de- 
termined from a numerical finite-element analysis. For 
this mode, the end faces at z = 0, I oscillate in oppo- 
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sition about the plane at z = 1/2 with frequency ujq. 
In thermal equilibrium at temperature T, Su z at the 
central points x, y — on the end faces has ampli- 
tude (SuDujq c± ksT /MqoUq, where Mo is the effective 
mass [12 j, which for the mode shown is Mq = 2.7 M. Sig- 
nificant to our current investigation is the axial strain per 
unit displacement, £ = e zz /u zi where £ ~ — 1600/m at 
z = for the parameters in Fig. [TJa) (i.e., expansion 
at z = with u z < generates strain e zz > 0). 

As depicted in Fig. [I] we next address the question 
of the phase shift for light reflected from the fluctu- 
ating surface of the substrate at z = 0. Following 
Levin [14], we introduce the displacement q(z) weighted 
by the normalized light intensity ip(r) over a plane at 
depth z > 0, namely q(z) = j dxdySu z (r, z)ifj(r), where 
ijj{r) = (2/ttWq) exp(— 2r 2 /wl) with r = ^ x 2 + y 2 . The 
incident field with vacuum wavevector k experiences a 
phase change SO = —2kq(z = 0) = —2kqo for the re- 
flected field due to the "piston" motion of the surface. 

In addition to SO, strain that accompanies surface mo- 
tion modifies the optical coating and thereby leads to a 
phase shift S/3, with S/3 expressed relative to the front 
surface of the coating. The overall phase shift for the 
reflected field is then S<$> = SO + 5(3, with \5(3\ < \S0\ for 
typical optical coatings. We now present new designs for 
the coating to achieve S<& — SO + S/3 ~ SO/Qq. 

The coating structure shown in Fig. [TJb) has an em- 
bedded resonant layer that functions as a FP cavity and 
gives rise to a rapid phase variation near the cavity res- 
onance [28]. Explicitly, the coating structure is specified 
as n (r]ir]2) l r]i(j x f] FP )(r]ir]2) p ~ l n s . Starting from the 
vacuum side no = 1, there are I double-layers (771 772 ) l 
with 771 5 2 = 7r/2 at the reference wavevector ko, followed 
by an n\ layer with 771 = tt/2 at ko, then the FP ri\ layer 
with single-pass phase shift j x rjpp at ko, followed by 
the terminating p — I double layers, and lastly the sub- 
strate with index n s . We restrict attention to the case 
of a "thin" coating for which the total coating thickness 
d <C a, wo and further assume that the axial strain in 
the coating e c zz is the same as that in the substrate e zz . 
These assumptions simplify the discussion of the essen- 
tial aspects of our scheme for noise compensation; see 
Refs. [23j [25] for detailed treatments of multilayer coat- 
ings. A discussion of thermal noise from the coating itself 
is deferred to our concluding remarks. 

Results for two particular coatings are given in Fig. [2] 
for wq <C a. For definiteness, we assume a coating stack 
made from layers of S1O2 with index n\ = 1.45 and 
TCI2O5 with ri2 = 2.03. Parts (a, b) are for displacement- 
induced strain £ = —1600/m, as appropriate to £ ( r ) i n 
Fig. [TJa). The coating structure is [p = 33, 1 = 8] with 
three curves drawn for mode orders j = 1,4,16 with 
t]fp = 7r. In (a), increased length for the resonant struc- 
ture leads to greater strain sensitivity, with \Sf3\ becom- 
ing larger relative to \S0\. For j = 16 there arise "magic" 
wave vectors k± for which 5&(k±) = 0, with then the 
piston phase shift SO from fluctuations of the surface dy- 
namically compensated by the strain-induced phase shift 
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FIG. 2: Phase shift S<f>(h) and transmission coefficient T(k) 
for two different coating designs illustrated in Fig. |TJb) and 
specified in the text, (a, b) The case £ = —1600/m for in- 
creasing phase jrjFP of the internal resonant structure, (a) 
Total phase 5<&(k) with "magic" wavevectors k± and (b) as- 
sociated transmission coefficient T(k) for j — 1,4,16 from 
widest to narrowest curves, (c, d) 5<&(k) and T(k) for the case 
£ = — 5000/m for two coatings with slightly different internal 
resonances 77^^ with distinct values fc^, k+ for nulling noise 
from mirrors B, A. In all cases, 5<f>(k)/\56(ko)\ is plotted for 
<7o > 0, with 50(h) given by the dashed lines in (a, c). 



Sf3 from the coating. At these "magic wavelengths," the 
phase of the reflected optical field becomes insensitive to 
the thermal motion of the surface of the substrate. 

Figures [2jc,d) investigate the possibility that thermal 
noise of individual mirrors of an interferometer could be 
measured in situ. The coating for each of two mirrors 
A, B is specified by [p = 33, 1 = 7, j = 8]. The thicknesses 
of the FP layers in the A, B coatings have been adjusted 
to give 77^ B) = (0.9995, 1.0005) x tt at k for the two 
traces shown. For operation at the lower (upper) value 
k^(k+) noise from mirror B(A) would be nulled while 
that from mirror A(B) would be ~ 80% of the full piston 
phase, so that the noise arising from mirrors A, B could 
be individually measured. For operation at the central 
value fc/fco = 1, thermal noise from both mirrors A, B 
would be suppressed. Panel (d) gives the transmission 
coefficient T(k) = 1 — R(k) for the two coatings, where 
T(kP,k+) is consistent with high-finesse measurements. 
Not shown in Fig. [2] are FP resonances away from k ~ fco, 
which can also be employed to suppress thermal noise. 

Because the surface strain |£| decreases with increasing 
size of the substrate, compensation of the piston phase 
SO with increasing a, wq requires greater departures from 
standard coating designs, material specifications, and 
fabrication procedures. Moreover, not all eigenmodes of 
oscillation have the proper parity (i.e., £ < 0) for com- 
pensation by way of the coating designs in Fig. [2] 

A second, more challenging example is thermal noise 
at frequencies uo well below the lowest resonance of the 
mirror substrate. In this quasi- static regime, we must un- 
derstand the correlation between thermal displacements 
at depth z and fluctuations of the surface arising from 
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many modes [12]. The foundation for our analysis is the 
FDT [Hl[9] as applied by Levin to this setting [14]. 

We consider a substrate in the form of an infinite half 
space with boundary at z = and extending to z > 0. 
Equation (9) in Mindlin [29] provides the Green's func- 
tion required to deduce the admittance for application 
of the FDT from Eq. (6.8) of Ref. [9]. Following the 
techniques in [16 , we find the spectral correlation for 
beam-averaged axial displacements q(z) at depths zi,z 2 
to be 



C( Zl ,z 2 ) 



(q(z 2 )q(z 1 )) UJ = (qoqo)ujN(z 2 ,z 1 ) 



(1) 



where q(zi) = q(zi,cu) is the Fourier transform of q(zi,t). 
(qoQo)cj — 2ksT(l — a 2 )(j) s /7r 3 / 2 woEu; is the standard re- 
sult for thermally driven surface displacements, with E, a 
the Young's modulus and Poisson ratio, respectively, and 
(j) s the (possibly frequency dependent) loss angle for the 
substrate. The function N(z 2 ,Z\) determines the me- 
chanical admittance and is given by 



N(z 2 ,z 1 ) 



where 



u) 



8^(1 - a) 2 J 



dke- k2 < /A f{z 2 ,z 1 -k), 



f(z 2 ,z 1 ]k) = e 



-k\Z-\ 



[3 — 4cr + fe|2r_|]- 



e~ hz+ [5 - 12a + 8a 2 + fe(3 - 4a)^ + + 2k 2 z 1 z 2 ] , (2) 

with z± = z\ zb z 2 and N(z, z) plotted in Fig. [3^ a). Note 
that iV(0,0) = 1, N(z 1 ,z 2 ) = N(z 2 ,z 1 ) [10], and that the 
admittance derived from N(z,0) agrees with that from 
previous work [14] [I6j [17] . 

Equation ([!]) enables us to determine spectral corre- 
lations between thermally driven axial displacements at 
depths z\ , z 2 characterized by [30] 



C(zi,z 2 ) 



(q(zi)q(z 2 )) 



(^O^i)) 1 / 2 ^)^))!^' 



(3) 



where G{z,z) = 1 and \C(z\,z 2 )\ < 1. From Fig. |3ja), 
we see that thermal fluctuations of q(z) correlate over 
length scales set by wo, albeit with decreasing amplitude 
(q(z)q(z)) oc N(z,z) away from the surface. 

Figure [3jb) investigates correlation between displace- 
ment q(zi) at z\ and axial strain at depth z 2 by way of 
the function Q(z 1 -,z 2 ,Az 2 ) = s c z oh (z u z 2 )/e z ot (z 2 , Az 2 ), 
with \Q\ < 1. Q(zi]Z 2 ,Az 2 ) expresses the ra- 
tio of coherent to total strain at z 2 within a small 
slice Az 2 , where e z oh (zi, z 2 ) = {q(z 1 )\q(z 2 + Az 2 ) - 
q(z 2 )]/ Az 2 ) / (q(zi) 2 ) 1 / 2 is the strain at z 2 correlated 
with the displacement q(z\) at z\, while s t z ot (z 2 ^ Az 2 ) = 
([q(z 2 + Az 2 ) — q^)} 2 ) 1 / 2 / Az 2 is the total strain at z 2 . 

In Fig. [3^b), the spatial scale for correlation of dis- 
placement and strain is again set by w$, but now with 
magnitude reduced by ^ Az 2 /wq. This scaling of Q 
can be understood from the fact that the rms strain 
diverges as 1/VV for thermal fluctuations in a vol- 
ume V [31], which motivates our use of finite differ- 
ences to characterize strain. Near the surface at z = 




FIG. 3: (a) Correlation of thermally driven axial dis- 
placements (7(21,22) and (b) displacement-strain correlation 
Q(zi]z 2 , A22) / \/ Az2 1 'w versus z 2 /w for zi/w = 0, 1, 3, 10, 
with AZ2/W0 — 10 -3 in (b). The dashed trace in (a) is 
N(z2,Z2). (c) Illustration of the mirror geometry discussed 
in the text, (d) Spectral density of phase fluctuations F(z 2 ) 
for the reflected field E r for the configuration in (c) for an 
Si02 substrate. From top to bottom, the curves are for 
a — 1.5, 0.3, 1.0, 0.7 with the blue trace for a — a m in overlay- 
ing the curve for a — 0.7. The dashed curve is from a simple 
model that incorporates incoherent contributions from trans- 
verse strains to 5(3. All curves are for a — 0.2. 



with Az/wq « 1, we find from Eq. that e t z ot (z ~ 
0, Az) ~ [(qoQo) / (Azwq)) 1 / 2 , where the relevant volume 
is V ~ ttWqAz. On the other hand, the coherent strain 

e coh = e c z oh (0,0) (qoQo) 1/2 /w , so that \e c z oh /s z ot \ 

scales as ^ Az/wq <C 1 near z = 0. 

Simply stated, thermally driven surface motion qo 
arises from strains within the substrate over distances 
z > wq. Hence, our previous strategy with a thin op- 
tical coating of thickness Az = d <C wo employed as a 
surface-strain sensor for a single eigenmode will compen- 
state only a small fraction of the phase shift SO = —2kqo 
from the surface motion in the quasi- static regime. In- 
stead, we must find a geometry for which 5(3 has contri- 
butions over z > wq sufficient to compensate SO. 

An initial mirror geometry that attempts to achieve 
such compensation is illustrated in Fig. 3^c). A thin di- 
electric coating with reflectivity Rq is deposited on the 
surface of the substrate at z = 0, and a second reflecting 
surface is embedded at z = z 2 <C nkwo 2 /2 with reflectiv- 
ity R 2 — > 1. In the static case, an incident field Ei lin- 
early polarized along x is reflected from this two-mirror 
geometry to give a field E r with E r /Ei = Ae lV ^\ where 
(p is the roundtrip, internal phase, and A = 1 for R 2 = 1. 
Thermally driven fluctuations lead to a phase shift 5& 
for E r , S<& = —2k [% + a(k) (q(z 2 ) — qo)] , where the pis- 
ton phase is SO = —2kqo and the interferometer phase 
is 5/3 = —2ka(k) (q(z 2 ) — %). a(k) expresses the sen- 
sitivity of the composite cavity Ro,R 2 to strain driven 
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changes in both length and index around ipo, with 

dT(tp ) 



a(k) = n s I 1 



dtp 



(4) 



where p\ 2 is an element of the strain-optic tensor pij for 
an assumed isotropic material. In Eq. [4j we neglect con- 
tributions from fluctuating strains transverse to the prop- 
agation direction z, which are considered below. 

The expressions for S& and a(k) lead to the spec- 
tral density ((S&) 2 ) of phase fluctuations for the re- 
flected field E r , namely ((S<$>) 2 ) = ((SO) 2 ) [(1 - a) 2 + 
2a(l-a)N(0,z 2 ) + a 2 N(z 2l z 2 )). Figure ^d) plots 
F(z 2 ) = ((S§) 2 ) I ((SO) 2 ) for several values of a. For 
fixed Rq (with R 2 = 1), a(k) varies periodically from 
minimum to maximum over the range Ak = it/n s z 2l 
with then a determined by the selection of k. F < 1 
represents phase noise reduced below that from ther- 
mal fluctuations in the piston phase ((SO) 2 ). Since 
S& = 2k [(1 — a)qo + aq(z 2 )] , a = 1 corrresponds to di- 
rect compensation of the piston phase SO oc % for the 
reflected field, albeit at the price of noise q(z 2 ) from fluc- 
tuations in the position of the reflecting surface at z 2 . 
More generally, the minimum Fmin for a given value of 
z 2 represents a compromise in noise from q and q(z 2 ) 
determined by the value a m i n . Fmin — 0.36 for z 2 ^> wq 
and a = 0.2 in Fig.^d). 

An important caveat related to Fig. |3|d) is that the 
full curves omit fluctuations in optical path arising from 
transverse strains e xx , e yy , which contribute by way of p^j 
to Sj3 and give rise to a scaling F ~ z 2 /wo for z 2 /wq ^> 1 . 
A full treatment of these effects is beyond the scope of 
our current analysis. Instead, the dashed curve in Fig. 
[3jd) is from a simple model based upon the FD theorem 
applied to e xx , e yy with loss angle (j) s and provides a rough 
estimate of their incoherent contribution to S(3. 

The conceptual design in Fig. 3^c) is likely far from 
optimal. Because the strain field associated with ip(r) 
at z = spreads transversely for z > 0, it is not 
well matched to our assumed optical profile with fixed 



w(z) = wo. Geometries with partially reflecting surfaces 
distributed along z might further reduce F • For finite 
thickness of the substrate, a treatment as in [16 is re- 
quired, with now the possibility of reflection from the 
rear surface [27] . More generally, coherent measurements 
over a range of k values could enhance sensitivity since 
S/3(k) can be tailored to be distinct from S0(k). 

Although our treatment has been exclusively for 
"Brownian" noise arising in the substrate, we suggest 
that our methods should be relevant to the suppression 
of thermal fluctuations from other sources within the sub- 
strate, such as thermoelastic-damping [I6j [THJ [19] . More- 
over, variations in the coating design from Fig. [2] could 
lead to schemes for suppression of thermal noise within 
the coating [23]. In contrast to the substrate for which 
strains at z <C wo have small correlation with §o i n the 
quasistatic regime, thermal noise from the coating leads 
to surface strains that are highly correlated with go- 

Certainly, important questions remain related to our 
proposals for noise compensation, including significant 
fabrication challenges, the impact of optical absorption 
within the coating and substrate, and the need for more 
complete theoretical analyses. We make no claim of a 
"magic bullet" for the elimination of thermally driven 
phase fluctuations in optical interferometry. Rather, our 
goal is to provide a perspective on thermal fluctuations 
which moves beyond a surface-centric view to consider 
the statistical character of the underlying stochastic dis- 
placements and strains that conspire to displace the sur- 
face and thereby to suggest new strategies for enhanced 
sensitivity and stability of optical interferometers. 
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